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Abstract 

The scattering is studied using moduli space metric for well-separated vortices of non- 
Abelian vortices in (2+l)-dimensional U(N) gauge theories with iV Higgs fields in the 
fundamental representation. Unlike vortices in the Abelian-Higgs model, dynamics of non- 
Abelian vortices has a lot of new features; The kinetic energy in real space can be transfered 
to that of internal orientational moduli and vice versa, the energy and charge transfer be- 
tween two vortices, the scattering angle of collisions with a fixed impact parameter depends 
on the internal orientations, and some resonances appear due to synchronization of the 
orientations. Scattering of dyonic non-Abelian vortices in a mass deformed theory is also 
studied. We find a bound state of two vortices moving along coils around a circle, like a 
loop of a phone code. 
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1 Introduction 



Topological solitons are localized finite energy solutions of classical equations of motion in field 
theories. Their stability is protected by topological winding numbers. Various topological solitons 
are known so far in field theories, for instance, instantons, magnetic monopoles, vortices and kinks 
(domain walls) [TJ. They are not only the classical solutions of the equations of motion but also 
play important roles in non-perturbative quantum effects through, for example, strong-weak 
dualities such as the electromagnetic dual in four dimensions. While topological solitons have 
been observed in condensed matter systems, in cosmology, solitons, especially cosmic strings, 
may be observed in a future by a direct detection of gravitational waves, the gravitational lens 
or the cosmic microwave background. 

In each dimension, we can consider dynamics of various topological solitons which behave as 
if they are particle-like objects ; kinks in d — 1 + 1, vortices in d — 2 + 1, magnetic monopoles 
in d = 3 + 1, and instantons in d — 4 + 1. In principle, the motion of solitons are determined 
by field equations which are usually highly nonlinear differential equations. Therefore a very 
complicated analysis is required to understand the soliton dynamics in general. For complicated 
physical systems, it is important to extract essential degrees of freedom by throwing away other 
unimportant degrees of freedom. This procedure is a challenging and interesting problem. Among 
various kinds of solitons, gauge theories often admit an important class of solitons, called local 
solitons; Yang-Mills instantons [2], 't Hooft-Polyakov monopoles [3], Abrikosov-Nielsen-Olesen 
(ANO) vortices jl] and CP 1 kinks (domain walls) |6]. A nice way of extracting such essential 
degrees of freedom has been established for a particularly important class of local solitons, called 
Bogomornyi-Prasad-Sommerfield (BPS) solitons [7|. BPS solitons saturate the lower energy 
bound and are the most stable among configurations with a fixed topology [7] . They further nat- 
urally appear in supersymmetric gauge theories, break/preserve a fraction of supersymmetry, and 
consequently are quantum mechanically stable under perturbative or non-perturbative quantum 
corrections [8]. Since static BPS solitons exert exactly no forces among them at any distance, 
multiple solitons can statically coexist at any position, which become parameters of multi-soliton 
solutions, namely moduli parameters associated with massless modes of solitons. Even though 
no static forces exist, solitons feel forces depending on their velocity. As long as solitons move 
slowly, it is sufficient to consider the motion of massless modes to describe the dynamics neglect- 
ing all the massive modes. This approximation is called the moduli space (geodesic, or Manton) 
approximation, which was proposed by Manton to discuss the dynamics of BPS monopoles [9j [1]. 
The moduli space dynamics is a good approximation when the kinetic energy of solitons is much 
smaller than any mass scales of the theory. According to the moduli space dynamics, solitons 
move along geodesies of the moduli space of the soliton, where the forces depending on velocities 
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of solitons are represented as geodesic forces. 

Unfortunately, it is not very easy to get the moduli space metric in general. Only in few 
cases, the metrics are known explicitly; the prime example is the Atiyah-Hitchin metric of k = 2 
BPS monopoles in the SU(2) gauge theory. For multiple monopoles, only asymptotic metrics 
were obtained when monopoles are well separated [H]. The moduli spaces of well-separated 
monopoles were also obtained in gauge theories with arbitrary gauge groups [12]. The moduli 
space approximation has been successfully applied to many other solitons such as ANO vortices 
in the BPS limit p31 UHl Ml El HE] , lumps pj], and BPS domain walls [6]. It has been applied 
even to BPS composite solitons of different kinds [201 121], such as domain wall networks (webs) 
[2"2"j [221 [21] and vortex-strings stretched between parallel domain walls [221 1231 121] • However 
only the asymptotic metric is explicitly known for well-separated ANO vortices [161 El 118] , 
because the ANO vortex equations are not integrableE These metrics have been used to discuss 
the scattering problems. In particular, scattering processes of solitons have attracted attention 
of many mathematicians and physists. It is well known that the magnetic monopoles scatter, 
surprisingly, with 90 degree, when they collide head-on [ID] . The same has been seen for the 
head-on collision of ANO vortices in the BPS limit [HI [15], and vortex-strings stretched between 
parallel domain walls |21j . 

In general the moduli space of BPS solitons is a surprisingly big space. Usually, the dimensions 
of the moduli space is proportional to the topological number (the number of the solitons) k. 
For example, the ANO vortices in the Abelian-Higgs model (Ginzburg-Landau model) at the 
critical coupling have 2k moduli parameters, which correspond to positions of the vortices [28] . 
This is much larger than the dimension of the symmetry group of the theory. Some topological 
solitons admit more moduli than the degrees of freedom of their positions. For example, the k 
BPS monopoles in the SU(2) gauge theory is known to have 4k = (3 + l)k degrees of freedom. 
As in the vortex case, 3k can be identified as the monopole positions. The remaining k degrees of 
freedom are U(l) phases of the internal space, which can be called as internal orientations. There 
are other examples of solitons which possess the U(l) orientational moduli; kinks in the CP 1 
model [5] and U(l) or U(N) gauge theories coupled to Higgs fields with non-degenerated masses 
[61 [2H]. Hence these solitons can be called Abelian solitons. On the other hand, non-Abelian 
moduli are associated with non-Abelian solitons; Yang-Mills instantons, non-Abelian monopoles 
|30j . non-Abelian vortices [3Tj, and non-Abelian kinks (in U(N) gauge theory coupled to Higgs 
fields with degenerated masses) [32]- The internal orientations can be regarded as Nambu- 
Goldstone zero modes corresponding to global symmetries or (global parts of) local symmetries, 

1 For ANO vortices in a hyperbolic plane with a particular curvature, the vortex equations become integrablc 
[26] and consequently the moduli space metric can be calculated [27] . 
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which are unbroken in the vacuum but are spontaneously broken in the presence of solitons. In 
general, the motion of internal U(l) orientations give preserved charges to the solitons, making 
them BPS dyonic solitons. 

Since the discovery of non-Abelian vortices [31] , much progress has been made in recent years 
[33l l23l [24]. Unlike the ANO vortices in the Abelian-Higgs model, the non-Abelian vortices have 
non-Abelian internal orientations and associated conserved charges; In the case of U(N) gauge 
theory with N flavors of Higgs fields in the fundamental representation, the internal orientation is 
the complex projective space CP^ -1 , which corresponds to Nambu-Goldstone modes associated 
with the SU(N)c + -p color-flavor locked global symmetry spontaneously broken in the presence of 
vortices. Because of non-Abelian internal orientations, we can expect that the dynamics of the 
non-Abelian vortices is much richer and more interesting compared to the ANO vortices, although 
the analysis gets much more complicated. The moduli space of multiple vortices with full moduli 
parameters was completely determined without metric by partially solving BPS vortex equations 

EU [Ml ES] ; The moduli space for k separated vortices is a /c-symmetric product 

M s k cp ~ (C x CP N ~ l ) k /S k cM k (1.1) 

of the single vortex moduli space [31] while the whole space M. k is regular. General formula for 
the moduli space metric and its Kahler potential were given in [36J . The metric of the moduli 
subspace for two coincident vortices [3HJ ESI EZ] , which is supplement to Ai s k e = 2 inside the whole 
space Ai k =2, was found, and it surprisingly shows that two non-Abelian vortices scatter with 
90 degree in head-on collision even though they have different internal orientations CP^ -1 as 
the initial conditions [3H] • Most recently, we have obtained the asymptotic metric on the moduli 
space M.\ ev of k well-separated non-Abelian vortices which is valid when the separation of vortices 
are much larger than the inverse Compton wave length of massive vector bosons, which is the 
length scale of the vortices [10] . 

In this paper, we study the dynamics of non-Abelian vortices by using the recently found 
asymptotic metric of non-Abelian vortices in the U (N) gauge theory with N Higgs scalar fields 
in the fundamental representation. In order to solve the dynamics, we will make use of the 
technique of the moduli space approximation. The asymptotic metric of the moduli space allows 
us to solve the dynamics of the well-separated and slowly moving non-Abelian vortices. We 
find that the dynamics of the non-Abelian vortex is quite different from that of the Abelian 
one. The major reason of the difference can be traced back to the conserved charges which 
are absent in the Abelian case. The vortices with the same charges repel while those with 
the opposite charges attract. Since the charges can change during the scattering process, non- 
Abelian vortices experience rich and subtle forces, which produce sometimes the counter-intuitive 
or unexpected dynamics. We find several new features of the dynamics of non-Abelian vortices: 
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i) the scattering angle depends on the internal orientation, especially parallel orientations give 
repulsion while anti-parallel orientations give attraction, ii) the energy of real and internal spaces 
can be transfered, iii) the energy and charge transfer between two vortices occur, and iv) some 
resonances appears due to synchronization of the orientations. 

We also study the dynamics of the dyonic non-Abelian vortices in the mass deformed theory 
|4Tj , with the method of the moduli space dynamics. A new feature in this case is that a potential 
term appears in the low energy effective action. Therefore, the solitons experience two forces: 
The one is the geometric force and the other is the potential force. In a special situation, the two 
dyonic vortices drift on a circular orbit. This fact strongly suggests that the non-Abelian dyonic 
vortices can have a bound state. 

The paper is organized as follows. In Section [2] we introduce the model which allows the 
BPS non-Abelian vortices and review briefly the low energy effective theory of the two non- 
Abelian vortices, and the moduli space dynamics. We also give the asymptotic metric on the 
moduli space and the geodesic equations for the well-separated vortices. In Sec. [31 we study the 
scattering of two non-Abelian vortices. After defining the Noether charges of vortices, we give 
typical examples of numerical solutions of geodesic motion on the moduli space in Subsec. 13.21 
We then analytically study properties of dynamics; the geodesic force in Subsec. 13.31 scattering 
of two vortices with a large impact parameter by free motion approximation in Subsec. 13.41 and 
dynamics with zero impact parameter in Subsec. 13.51 In Sec. HJ we consider a mass deformation 
of the theory and investigate the dynamics of the dyonic non-Abelian vortices. We find a bound 
state of two dyonic vortices moving along coils around a circle. Sec. |5] is devoted for conclusion 
and discussion. In Appendix A, the effective action of vortices is written in the U(2) case in 
terms of unit three-component vectors. 

2 Asymptotic metric for non-Abelian vortices 
2.1 Lagrangian and BPS equations 

We consider a U(N) gauge theory in (2 + l)-dimensional spacetime with gauge fields for U (1) 
and W£ (a = 1, . . . , N 2 - 1) for SU(N) C , which couple to N Higgs fields H A (A—1,...,N) in 
the fundamental representation of the SU(N) C gauge group. The Lagrangian is given by 



C 



1 

4^ 



(2.1) 



V 




(2.2) 
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where £ is the Fayet-Iliopoulos parameter, e and g are gauge coupling constants for U(l) and 
SU(N) C , respectively. The overall scalar coupling constants in the potential V are chosen to 
be equal to the square of the gauge coupling constants, so that the model admits the BPS 
non-Abelian vortices. Thus the model has the three coupling constants e,g,£. In the three 
dimensional spacetime, all of the mass dimensions of e 2 , g 2 , £ are unity. Our convention is T>^H A = 
(d^ + iWfJP + iW^t a )H A and j '^t ' + F^ u t a = — z[£> M , V v \. The matrices t° and t a are the generators 
of U(l) and SU(N)c, normalized as 

t° = ^=l*r, Tr(t^) = V (2.3) 

As is well known, the Lagrangian Eq. (12. ip can be embedded into a supersymmetric theory with 
eight supercharges. The Higgs fields can also be expressed as an iV-by-iV matrix on which the 
SU(N)c gauge transformations act from the left and the SU(N)p flavor symmetry acts from the 
right 

H -> UcHU)?, U c G SU(N) C , U F e SU(N) F . (2.4) 

The vacuum of this model (V = 0) is in an SU(N)c+f color-flavor locking phase (U F = Uc), 
where the vacuum expectation values (VEVs) of the Higgs fields are 

1/2 

In this vacuum, we have a mass gap with the mass m e for singlet and m g for adjoint representa- 
tions of SU(N) c +f 

m e = e-y/c, vn g = gVc. (2.6) 
The energy density for a static configuration can be rewritten as 



H=yfcl N , c=(4) £■ (2-5) 



E = — Tr 

* 2e 2 



f 12 t°-e 2 (H{t°H A -0\ +^Tr F&P - g\H\t a H 



2 1 



ia ±a 2 / rrt j.a ttA\ 
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+A\V 2 H A \ 2 - £Tr[/ 12 t u ] - ie^d t (H\V 3 H A ). (2.7) 

For configurations with vorticity k (vortex number), the energy of is bounded from below by the 
following BPS bound 

E > kM v = -^ J d 2 xTr[f 12 t } = 2nck, fceZ, (2.8) 

where we have assume that the last term in Eq. (12.71) vanishes at infinity. This bound is saturated 
if the following BPS equations are satisfied: 

V- Z H = 0, \f 12 t° + ^-F? 2 t a = HH^ - cl N , (2.9) 

e z g z 



where z = x 1 + ix 2 is a complex coordinate. One can easily verify that all the solutions of the BPS 
equations solve the original equations of motion of the Lagrangian in Eq. (12. ip . The integration 
constants (the moduli parameters or the collective coordinates) contained in the solutions of the 
BPS equations parameterize the set of configurations with degenerate energy, that is, the moduli 
space of BPS vortices Aik- There are N complex moduli parameters for each vortex: one of 
N is position Zi and the rest N — 1 are internal orientations 0j (I = 1, 2, • • ■ ,k). Since no net 
forces are exerted among static vortices, each vortex has the position moduli zj as its degree of 
freedom. The internal orientation /3j is associated with the SU(N)c+f color-flavor symmetry, 
broken by each vortex down to SU(N — 1) x £7(1). The Nambu-Goldstone zero modes localize 
on each vortex and the corresponding moduli (3j parameterize the coset 

SU(N) ^ nJV _! 



CP N ~ L . (2.10) 



SU(N-l) x U(l) 

In the following, the (N — l)-dimensional vector 0i denotes the inhomogeneous coordinates of 
CjpTV-i for the internal orientation of 7-th vortex. The moduli space M.^ is a /ciV-dimensional 
Kahler manifold parameterized by the holomorphic coordinates zi and Pi and has the SU(N) 
isometry acting on f3j, which descends from the SU(N)c+f global symmetry in the vacuum. 



2.2 Asymptotic metric for non-Abelian vortices 

The low-energy dynamics of the vortex system can be described by an effective Lagrangian in 
which the moduli parameters are promoted to dynamical variables. The effective Lagrangian for 
these moduli parameters is given in terms of the Kahler metric on the moduli space 



L = gqW ', = 1, • ■ • , dim c M k = kN), 

where <p % are the holomorphic coordinates on the moduli space: 

{^} = {z I J I }. 



(2.11) 



(2.12) 



The metric gq of the moduli space consists of the free part and the interaction part which are 
given in terms of the corresponding Kahler potentials 



O 2 
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Khee + 



o 2 



int- 



d(f) l d(f>i d(f) l d(f)^ 
The free part describes the dynamics of completely isolated vorticejf] 



k r 



K 



free 



E 

i=i 



-M yZl 2 + — log 1 
2 g 2 



(2.13) 



(2.14) 



2 The Kahler class iw/g 2 can be determined [42] from the fact that sigma model instantons inside a vortex 
worldsheet are Yang-Mills instantons from the bulk point of view fJ3] . 
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c 







1.1363(7) 


0.25 




1.1853(1) 


0.5 




1.3090(5) 


0.75 


2.1955(9) 


1.48517(9) 


1 


1.7078(6) 


1.7078(6) 


1.5 


1.4714(7) 


2.3031(0) 


2 


1.4036(9) 


3.14(5) 


2.5 


1.3746(1) 


4.31(8) 


3 


1.3594(3) 


5.9(5) 


oo 


1.3266(7) 





Table 1: Numerical data for k = 1 U(2) vortex. The constants c e and c g are not well-defined for 
m e < 2m g and m g — > oo, respectively |44j . 



where M v = 2nc is the tension of the vortex and 4-7T / ' g 2 corresponds to the radius of CP^" 1 . On 
the other hand, the interaction part of the Kahler potential describes the leading interactions 
between well-separated vortices 



-2ttN 



KJ 



c 2 c 2 
-^K (m e \zij\) + -^QijKofaglzu]) 
^ 9 



(2.15) 



where K stands for the modified Bessel function of the second kind. The interaction term K^ 1 '^ 
between 7-th and J-th vortices is a function of the relative distance and an SU(N) invariant 
quantity 9/j defined, respectively, by 



\Zij\ = \zi - Zj\ 



|1 + /3J-/3j| 



Qu = N 1 Z HI H 1. (2.16) 

(l + l/^Xi + l/^l 2 ) 

The origin of the modified Bessel function K appearing in the interaction term can be traced 
back to the asymptotic tail of a profile function of the vortex. Since the leading term in the 
modified Bessel function Ko(m\zu\) is of order e _m ' z// ', the interactions exponentially vanish for 
large \zu\. Note that the mass scales m e and m g can be interpreted as inverse widths of "Abelian 
core" and "non-Abelian core" of a vortex, respectively. The strength of the asymptotic coupling 
is controlled by the constants c e and c g , which depend on the ratio m g /m e and N [S]. The 
numerical values of c e and c g for U(2) vortices are given in Table [TJ 

For later convenience, let us factor out 47r/<? 2 from the Kahler potential and rewrite it in 
terms of three coupling constants g, m e and m g instead of e, g, c 

k 



Kfree - — T 

q z 

y 7=1 



-m^ 7 | 2 + log(l + |/3 7 | 



Air 



rat 



E 

KJ 



N 2 

m n 

2 9 



mi 



; Ko(m e \zij\) + -^QijK^mglzu]) 



(2.17) 
(2.18) 
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2.3 Equation of motion along geodesies 



The overall coefficient 1/g 2 of the effective Lagrangian plays a role of the Planck constant (loop 
counting parameter) for fixed mass scales m g and m e . Therefore, the quantum effects are negli- 
gible and the classical analysis of the effective Lagrangian is valid for the energy scale E much 
larger than g 2 

g 2 <i?< m g , m e . (2.19) 

Note that the energy scale E should be much smaller than mass gaps to justify the use of the low- 
energy effective Lagrangian. In the classical analysis, the dynamics are completely independent 
of the gauge coupling g. 

The equations of motion for the moduli parameters <fi l take the following form of the geodesic 
equation 

l + l>'o A ' = 0, T) k ee g%g kI . (2.20) 

Since the interactions are sufficiently small (K^' J ^ « e - m \ z u\^ f or well-separated vortices, the 
Christoffel symbol Y l - k can be approximated as 

T) k « t) k + fV 3 d k d T K mt , (2.21) 

where Y l - k1 gq and V, are the free part of the Christoffel symbol, metric and covariant derivative 
respectively. Then, the equations of motion for zj can be written as 

zj = AETF^^' ( 2 - 22 ) 

m 9 1^1 " ZlJ 

where we have defined a differential operator 5u by 

d 2 d ( ? d ^ d \ d z d 



S„ s itj — + 2iu — ^p,. — + pj. —J + W , . w *. Wj . (2.23) 
The equations of motion for the orientations are 



1 + 1/3 



J+i L 



-^ + 8 I 8l- —^7 8u ^—Rj . — 

d(3\ fIfl dfi 1 + N 2 dPj 



K {I ' J \ (2.24) 



N-l 



where Vt/3/ is the free part of the covariant derivative along the trajectory on CP 

Vjj = & ■ 0i) (2-25) 

1 + \Pi\ 2 

We refer the right-hand sides of the equations of motion (12.221) and (I2.24p as forces, more precisely, 
geodesic forces. In general, they are proportional to squares of velocities oc d t 4> l d t 4^ ■ 



8 



3 Scattering of U(2) non-Abelian vortices 



In this section, we discuss the asymptotic scattering of the non-Abelian vortices in the simplest 
example of the U(2) case which shows essential differences between the Abelian and the non- 
Abelian vortices. In this case, the orientational moduli space becomes a sphere S 2 ~ CP 1 . The 
inhomogeneous coordinate fli is given by the stereographic projection from the north pole of S 2 
and related to the standard spherical coordinates as 



One can also make use of the following three-dimensional unit vector as coordinates of S 2 



Since 5*0(3) ~ SU{2)/Tj2 symmetry is manifest on Hi, the dynamics of the orientational modes 
can be better understood in terms of Hi. The Lagrangian written in terms of Hi is given in 
Appendix [A] 

3.1 Noether charge 

One of the sharp contrast of the dynamics of non-Abelian vortex to that of Abelian vortex 
comes from the internal orientations. As mentioned in the introduction, such internal degrees 
of freedom leads to conserved Noether charges. So let us begin with describing the conserved 
charges of non-Abelian vortices. 

Corresponding to the SU(2) global symmetry, we have one set of conserved charges Q. In 
the original theory, they are given by 




(3.1) 



Hi = (sin 9i cos <pi , sin 9i sin <pi , cos 6i). 



(3.2) 




(3.3) 



where a are the Pauli matrices. Although Q = for static non-Abelian vortices, these charges 
arise from the motion of the orientational modes. 



In the effective action of vortices, the charges of vortices are given by 




where £ 7 are the holomorphic Killing vectors associated with SU(2) rotation on the moduli space. 
Explicit form of the SU(2) triplet vector is given by 
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Fig. 1: The free motion of the internal orientation on CP 1 ~ S 2 . 

These charges have contributions from the free and interaction parts of the effective Lagrangian. If 
all the vortices are isolated, contributions from the free part of the Lagrangian for each individual 
vortex are separately conserved. In terms of Hi, we define the charge of the I-th vortex b)j^ 

Qi = -^ni x ni. (3.6) 

r 

For the isolated vortices, the dynamics of the orientation is described by free equation of motion, 
which is nothing but the conservation law of the charge of each individual vortex 

dt g 2 

This equation is equivalent to the geodesic equation on a sphere, so that the trajectory of the 
orientation is a great circle, as shown in Fig. [TJ 

ni = nio cos(c<;/i) + qi x uiq sin(a;/t), (3.8) 

where n/o and qj are unit constant vectors satisfying n/o ■ qi = 0. The charge of this solution is 
given by 

Q, = (3.9) 
T 

Note that we can always set uji > without loss of generality. 



iQi = —ni xHj = 0. (3.7) 



3.2 Features of the scattering of non-Abelian vortices 

As stressed in the previous subsection, the internal orientations and the conserved Noether 
charges make the dynamics of non-Abelian vortices quite different from that of the well-known 

3 The normalization of Qi is chosen so that they obey the half-integer quantization condition in the quantum 
theory. Although the coupling constant g appears in the classical equations of motion because of this normaliza- 
tion, it can be absorbed by rescaling Q] appropriately. 
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ANO vortices. Since the internal orientations are continuous parameters, there exist continu- 
ously different initial conditions, which can make the dynamics quite complicated. In order to 
see how much the dynamics of non-Abelian vortices differs from that of the ANO vortices, we 
show several numerical solutions for the geodesic equations (I2.22p and (I2.24p . 

The relative internal orientations between the first and second vortices are denoted as AO = 
9i — 9 2 , Acp — cpx — cp 2 - We specify initial conditions at sufficiently past. We can always use the 
SU (2) symmetry to set the initial values at sufficiently past as 

<^10 = -<P20, 010 = 020 = 0. (3.10) 

This initial condition corresponds to a pair of vortices whose orientations are rotating around 
great circles at the longitudes <^ 10 and </?2o — ~V?io- Therefore, we can choose five initial conditions 
for the orientations: &i (I = 1, 2) and Aip. 

Scattering of vortices without initial Q-charges 

As the first example of the scattering of two non-Abelian vortices, we examine the case of 
m e = m g and choose the following initial conditions for the relative orientation Atpo and for the 
velocities of the orientations 8io'- 

A<p = 0, 9io = 0, (1 = 1,2). (3.11) 

In this setting, the vortices have no initial charges Qjq = 0. Because of these initial conditions, 
the variation of the internal orientation is very small before the vortices passes through the 
interaction region. 

In Fig. (a), we show the scattering orbits for the initial relative orientation A# = tt/2 
by changing the impact parameter a as a = 2.0, 2.5, 3.0, • • • . One can clearly see that the 
moving vortices feel the repulsive geodesic force between each other. At a glance, this result 
appears to be very similar to the scattering of the ANO vortices, in which case the scattering is 
repulsive and the orbits are uniquely determined if the initial velocity and the impact parameter 
are fixed. However, it is actually different. See Fig. |2]-(b) where the impact parameter is fixed 
at a = 3 and the initial relative orientation is varied as A9q = 0, |, |, ^,7r. When A9 = 
(A# = j, |...), they scatter (almost) in the same way as the Abelian vortices. On the other 
hand, non-Abelian vortices just pass through each other without feeling any interactions for 
A#o = vr (a pair of vortices with anti-parallel orientations). This behavior is markedly different 
from the scattering of ANO vortices. We thus have found that the scattering angle is sensitive 
to the relative orientation. 

The dependence on the relative orientation can be roughly understood from the fact that 
the interaction part of the Kahler potential K int given in Eq. ( I2.18P is proportional to 1 + ©i 2 
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(a) a = 2.0, 2.5,---, A9 = f (b) a = 3.0, A0 O = 0, \ , \ , f , tt 

Fig. 2: Scattering orbits of non-Abelian vortices on the z-plane ((x, y)-plane) in the case of 
m e = m g = 1 (c e = c g = 1.708), for (a) various values of the impact parameter a with the 
fixed initial relative orientation A6 = f , and (b) various values of the initial relative orientation 
A#o with the fixed impact parameter a = 3.0. In the both cases the initial velocities of the 
orientations are zero (Qj = 0). No interaction exists for A# = tt (anti-parallel orientations). 
Otherwise, the interaction is always repulsive with reaching the maximum at A6q = (parallel 
orientations), where the interaction reduces to that of the ANO vortices. 

for m g = m e (c g = c e ). Since the relative orientation is almost unchanged from the initial 
condition until the vortices leave the interaction region, it follows that the scattering angle is 
also proportional to 1 + ©12. For example, the scattering angle is maximized (vanishes) for 
A# = 0, (AO = it) at which 1 + 612 = 2, (1 + 12 = 0). We emphasize again that, in contrast 
to the non-Abelian vortex scattering, the scattering of the ANO vortices is uniquely determined 
with the initial velocity and the impact parameter fixed. 

In this example, we have taken the initial condition without initial Q-charges. The numerical 
calculations show that the vortices are charged Q\ = —Q2 7^ after the scattering, even if 
orientational moduli is initially static. We will see this phenomenon via an analytic discussion 
in section 13.41 

Scattering of vortices with maximized non-Abelian effect 

As the second example let us consider the case where non-Abelian effects are maximized. If 
the Abelian vector boson mass is sent to infinity^ m e — > 00, the Abelian part of the interaction 

4 Our original model reduces in the limit m e — > 00 to a CP 1 nonlinear sigma model whose SU(2) isometry is 
gauged. 
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(a) a = 2.0, 2.5,---, A9 = tt (b) a = 3.0, A6 = 0, \ , f , f , tt 

Fig. 3: Scattering orbits of non-Abelian vortices in the z-plane ((x, y)-plane) for m e = oo, m 9 = 1 
(c g = 1.136), with (a) various values of the impact parameter a and the fixed initial relative 
orientation A^o = tt, and (b) various values of the initial relative orientation A6q and the fixed 
impact parameter a = 3.0. In the both cases the other initial conditions are: Atp = and 8j = 0. 

Kahler potential K[ nt , given by the first term in Eq. f 1 2 . 1 8 [) . is highly suppressed and only the 
non-Abelian part, the second term, survives. The numerical results are shown in Fig. |3j Except 
for the Abelian mass m e , the other parameters and the initial conditions are chosen to be the 
same as those for Fig. [2j Fig. |3]-(a) shows the scattering orbit for A9 = tt where the vortices 
clearly attract each other. This attractive force is a characteristic property of the non-Abelian 
case, which has not been seen in the Abelian case. As before, the scattering angle is affected by 
the initial relative angle of the internal orientations, as shown in Fig. [3]-(b). In the present case, 
the interaction is proportional to O12 = {1,0,-1} for A6 = {0,^,^}, respectively. Thus, the 
scattering angle changes its sign at A# = f • For A8 Q < |, the vortices repel as in the Abelian 
case, while they attract for A^q > |. The parallel orientation (A# = 0) gives the maximal 
repulsion and the anti-parallel orientation (A# = tt) gives the maximal attraction. 

Scattering of vortices with non-zero Q-charges 

The last interesting example is the scattering of non-Abelian vortices which have non-zero 
initial charges Qi Since a repulsive (attractive) force works between the same (opposite) charges, 
the scattering is very different from that with vanishing charges. The results are shown in Fig. m 
where 7/ = 6jo/(m g v) is defined by the initial relative velocity v and the angular velocities 9iq. 
Fig. 0]- (a) shows the scattering orbits of vortices with the same initial Q-charges. Because of the 
repulsive force, the vortices recoil even for a = 0. On the other hand, the vortices with opposite 
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(a) a = 0.0, 0.5, • ■ ■ , A(p = n (b) a = 4.0, 4.5, • ■ ■ , Ay? = 

Fig. 4: Scattering orbits of non-Abelian vortices with non-vanishing conserved charges for 7! = 
~l2 = f, rn e — m g — 1 (c e = c g = 1.708), with various values of the impact parameter a and of 
the initial relative orientation Atp . The other initial conditions are : #10 = #20 — f • 



charges feel attractive force as shown in Fig.0]-(b). We can see that their orbits are slightly 
wavy. This is due to oscillations of the forces caused by the rotations of the orientations. In the 
next section, we will discuss an "average" over the rapid motions of the orientations to find out 
effective forces between non-Abelian vortices with Q-charges. 



3.3 Forces between non-Abelian vortices 



Now let us examine the geodesic forces between two vortices induced by motions of the moduli 
parameters. To this end, let us first rewrite the equations of motion Eqs. (I2.22p and f)2.24p in 
terms of ni asf| 

„ <■» 

Z\2 



c c 

—K^rrielzul) + —(ni ■ n 2 )Ki{mg\zi 2 \) 
m e m g 



^12 / • n2 

{m g z l2 ) 



\ Z 12\ 



2c 2 q 1c\ Zl2 
H -K (m g \zi 2 \) (ni • a 2 + n 2 ■ a x ) m g z 12 -Ki(m g \zi 2 \)di ■ a 2 —— v (3.12) 



Z 12 



5 Note that the equation of motion n x n = A is equivalent ton = ixn - n\n\ 2 for vectors n and A such that 
|?i| 2 = 1, n- A = 0. 
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rii x ni 



Re 



-c 2 g K 2 (m g \z 12 \) n x x (n 2 - ffli x n 2 



+ 2c Ki{m g \zi2\) iol\ \ni ■ n 2 prrr- , \rn g z 12 



- c g [ K (m g \zi 2 \)iai(n 2 -ai + 2n 1 -a 2 
where a.j (7 = 1,2) is a complex three- vector defined by 




m g zi 2 



(3.13) 



ati = ni — in i x nj. 



(3.14) 



Note that the equation of motion for n 2 can be obtained by exchanging n\ and n 2 in Eq. (I3.13p . 

The geodesic force given in the right-hand side of Eq. (13.1 2p shows that the motion in orien- 
tational moduli space generally induces a force between the vortices even if there is no spatial 
motion initially. Similarly, Eq. (13. 13[) shows that a spatial motion induces a motion in orienta- 
tional moduli space, even if orientational moduli is initially static. This implies that the kinetic 
energy in spatial motion and in orientational moduli transmute each other. 



3.3.1 Average over rapid motion of internal orientations 

As we have seen in the numerical example Fig.|4]-(b), the sign and magnitude of forces between 
them oscillate due to the rotations of nj. Since motion in orientational moduli is (almost) 
periodic, its physical effect is best seen by averaging over the periods of two orientational moduli 
individually. Denoting this time average for the high-frequency modes by ( ) , we can express the 
averaged equation of motion for the relative position z 12 as 

(2ia) = ( ~r*7 \clm e K x {m e \z 12 \) ( ^fX - ^-K x {m g \z 12 \) Q x ■ Q 2 ) V (3.15) 
\ 1*12 1 \ V \ Z 12\ J 2n 2 m g J J 

This approximation is valid only when the motion in space is slow compared to the velocities of 
the orientations, 

cux, uj 2 3> Tn g \z\2\ and u)\ ^ u 2 . (3.16) 

. The first term depends only on the velocities of the vortices just in the same way as the force 
between Abelian vortices. The second term is the dominant force induced by the motion of 
orientational moduli. Note that the order of the first term is m e \zi 2 \ 2 , while that of the second 
term is ujiu 2 /m g . Thus, in the case of m e > m g , the first term is negligible compared to the 
second term under the condition (I3.16p . On the other hand, if m e < m g , the first term can be 
dominant since K (m e \zi 2 \) 3> Ko(m g \zi 2 \) asymptotically. 
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It has been shown that the interaction of non-Abelian vortices are well described by regarding 
them as point-like sources of the Higgs fields and massive vector fields [3D]. The motion of the 
orientation of the J-th vortex induces the following (color) electric charge distribution^ and the 
massive vector field 

fo = ^ Qi ■ B d z d- z 5 2 (z - Zl ), Wi = ^Qj-B K (m g \z - z T \). (3.17) 



Therefore, the electrostatic potential between J-th and J-th vortices is given by 



j d 2 xTi 


— j J W J 
g2 Jo vv o 







2 2 

^-Qi-Qj K (m g \ Zl j\). (3.18) 

The second term in the averaged equation of motion Eq. (13.151) can be attributed to this "elec- 
trostatic interaction". This force is repulsive (attractive) if the inner product of the charges of 
two vortices is positive (negative). Namely the vortices with aligned charges repel each other, 
whereas those with disaligned charges attract each other. These phenomena can be seen in the 
numerical calculations in Fig.|U 

Let us next investigate the averaged equations of motion for the orientations. By using the 
charge vectors Qj defined in Eq. (13.61) . the equation of motion of ffi Eq. (I3.13P reduces to 

^ = ^ (Ko(m g \z 12 \) (Q 2 x &)} , (3.19) 

where we have used (nj ® nj) = \(l — qi ® qf}. The equation for Q 2 can be obtained just by 
exchanging the two vortices (1 -H- 2). These equations describe precessions of the orientations, 
that is, motions along great circles with slowly moving axes. We can see from Eq. (I3.19P that 
the sum Qtot = (Qi + Q2J an d the inner product (Qi ■ Q 2 j are conserved. By solving averaged 
equation of motion, we find that each charge Qi slowly rotates around Qtot as (see Fig. [5]) 

Qi) ^ a Qtot + Qosc cos 6(t) + Qto * * Qosc sin 9(t) , (3.20) 

I Qtot I 

Q 2 ) ^ 6 Qtot ~ Qosccosfl(t) - gt ^ X g ° sc sing(t), (3.21) 

I Qtot I 

where a + b = 1, Q osc ■ Q tot = and 

= ^\Qtot\ f dt K (m g \z 12 \). (3.22) 



We can see that the angular frequency of the rotation is quite small 9 0J\ y2 . 



6 The charges Qi themselves are not the electric charges but the charges of SU (2) global symmetry. The total 
electric charges are always zero even for vortices with rotating orientations. 
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. Qtot 




— * 

Fig. 5: An example of motions of the charge vectors Qi for uj ^> va g \z\2\- 

In contrast to the Abelian vortex, the non-Abelian vortex has the characteristic feature that 
the interactions depend on the relative internal orientations. Since the inner product Qi ■ Q2 
is almost conserved during the scattering process, the interaction can be, for instance, kept to 
be attractive by choosing near anti-parallel orientations. In such a situation, one can imagine 
a possibility for two non-Abelian vortices to be bound together. However, it is not the case. It 
is because the balanced point for the exponentially decaying attractive force (~ e m a\ z ^ with 
the centrifugal force (~ l/|z 12 | 2 ) is unstable. Of course, this conclusion is valid only for the 
well-separated vortices at distances \zi 2 \ > l/m g . Although one may expect that bound states 
can exist when \z\2\ < l/ m 35 such a region is out of range of validity of our approximation. In 
Sec. HI we will see that a bound state of two (dyonic) vortices exists in a mass deformed theory 
where a certain extra force, which is not a geodesic force, is induced by a mass term. 



3.3.2 Limit of slowly moving internal orientations 

Let us next consider the opposite limit of Eq. (I3.16p . 

wi <^ m g \zi 2 \- (3.23) 

Now the motion of the internal orientations is very slow compared to the spatial motion, so that 
we can safely neglect Hi in the equations. Thus the first lines of Eq. (I3.12p and Eq. (13 . 1 3 [) give 
dominant contributions. Depending on the mass scales m g and m e , the dominant contributions 
in Eq. (13. 121) are given by 

--^-K^m^z^j— ;(m e z 12 ) 2 for m e < m g 



Z12 



m e \z 12 \ 

„2 



,,2 



*12 1) 


Z12 1 


— r( m 


^12 1 








{m g z 12 ) 


1^12 1 
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for m e = m g (3.24) 
for m P > m„ 



Since the interaction for m e < m g is equivalent to that of the Abelian case, the vortices receive 
repulsive force in the scattering process. For m e > m g , the forces depend on Hi ■ n 2 . When 
Hi ■ n 2 > (Hi ■ n 2 < 0), they repel (attract) each other and the interaction accidentally vanishes 
at Hi • H 2 = 0. Although the interaction depends on Hi ■ n 2 also in the case of m e = m g , it is 
qualitatively similar to the Abelian case, that is, the vortices receive the repulsive force except 
for the case Hi ■ n 2 = — 1. This explains the behaviors of the scattering given in Figs. [2] and [31 

3.3.3 Internal orientations at a head-on collision 

As a special case, let us consider a pair of vortices which are going to collide head-on with each 
other 

Z\ 2 — — w-— ni,n 2 ~0, v > 0. (3.25) 

1^12 1 

From the equation of motion for the orientation Eq. f)3.13p . we find that 

Hi x Hi ~ c 2 g m 2 g v 2 K 2 (m g \z 12 \) n Y x n 2 . (3.26) 

This equation implies that the orientations of the vortices tend to align before a head-on collision. 
This result for well-separated vortices naturally extends a similar result obtained by an analysis 
around a vicinity of coincident vortices |39j . 

3.4 Large impact parameter 

If two vortices are completely separated, interactions between them can be neglected, and the 
solution of the equations of motion is given by the free motion (moving parallel to the x-axes) 

z l2 = vt + ia, n 7 = c^' 1 + c/e _iWl *. (3.27) 

Here we have rewritten the solution Eq. (13.81) in terms of a complex vector cj which is related to 
Hio and qi as 

ci = - (H I0 - iqi x n IQ ) ^ H I0 = c 7 + c/, <f/ = -2ic 7 x c/. (3.28) 

Let us here consider the case of large impact parameter a. If two vortices are far apart, their 
interactions are exponentially suppressed and the deviation from the free motion is small during 
the scattering process. Therefore we can safely evaluate their interactions by approximating the 
right-hand side of the equations of motion in Eqs. f)3.12p and (I3.13P by inserting the free motion 
in Eq. (153711 . 
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3.4.1 Scattering angle 



Substituting the free motion Eq. (I3.27P into the right-hand side of Eq. (I3.12p . we obtain the 
approximated equation of motion for the relative position z\ 2 in which all the forces are known 
functions of time. Then, we can evaluate the total change in the relative velocity by integrating 
the forces from t = — oo to t = oo 



Aii 2 



Z\2 dt. 



(3.29) 



This can be carried out explicitly by using the following formulas of Fourier transformations 

71 



dtK (my/(vt) 2 + a 2 )e 



ILUt 



a/ (mv) 2 + a; 2 



dt / Vt±m K^my/ivty + a^ = — 
^(vt) 2 + a 2 mv 



a/ (mv) 2 + ui 2 

Eventually, we obtain the following deviation of the relative velocity 
Ai 12 1 



•fV(™>) a -K» a 



± 1 



(3.30) 



■§V( nM ') 3+fc ' a .(3.31) 



2niv 



-de mea + d 7 i7 2 (gi-g 2 )e- 



(3.32) 



-c„72 



1 + it 02 + ill Qi x q 2 • nio e 



ji\/l+7i 



_C 9^ 1 1 + 72 Ql ~ il2 Ql X 92 J ' ™ 20 e ' 

71 + 72 



+c 2 (l + 2 7l72 ) Re(ci • c 2 ) 



a/1 + ( 7 i + 72 ) 5 



:im(ci • c 2 )j e- m » a Vi+W 



+4(1 - 2 7l72 ) ( Re(c 1 • c 2 *) + 71 72 ilmfo • c 2 *) ) e^V^^ 3 ^ 2 

a/I + ( 7l - 72 ) 2 



Note that the frequencies uj appear only through the following ratio 



li 



TflgV 



(3.33) 



This follows from the invariance of the ratio Az^/v under the time rescaling t — > At, which arises 
from the property of the geodesic equation (I2.20p itself. This ratio 7/ measures how much the 
vectors nj rotate in the typical time scale of the interaction At m l/m g v. For generic values of 
the ratios 7/ and the vectors qj, only the first line in Eq. ( 13 .32 j) is dominant, 



7iivc 2 e mea 



(3.34) 



27nwc 2 7i7 2 (<fi • q 2 ) 

In the case of m e < m g , the leading contribution comes from the first term in the averaged 
equation of motion Eq. (I3.15p . while the second term is dominant for m g < m e . 



19 



The free-motion approximation is valid for orbits with small scattering angles. This means 
that the kinetic energy of the spatial motion should be sufficiently larger than that of the internal 
orientations. If vortices with non-zero charges have small relative velocities, the interaction 
Eq. (I3.18P lasts for a long time interval At rs l/(m g v) and the deviation of the trajectory from 
the free motion becomes large. Therefore 71 and 72 should be generically small so that 

7i7 2 e- m9a < 1. (3.35) 

In Fig. |6]-(a), we compare the scattering angles obtained by numerical calculations and those 




X 



(a) A% (b) Orbits in the z-plane 

Fig. 6: (a) Scattering angles as a function of 72, and (b) the corresponding orbits of the vortices in 
the z-plane ((x, y)-plane) for m e = m g = 1, a = 7, q± = q% = (1,0,0), nio = (0, 1,0), = (0,0, 1) and 
7i = 3. 

with the free-motion approximation. One finds that the free-motion approximation is indeed 
valid only for the scattering with small scattering angles. We also show the numerical results 
of the scattering with various initial conditions in Fig. [6]-(b). Among them, there are orbits 
whose scattering angles exceed ir/2. For such collisions, the free-motion approximation cannot 
be applied. 

The terms with 7/ in the exponents (terms other than the first line) in Eq. ( 13.321) are contribu- 
tions from the oscillating forces and become smaller for larger values of 7/. Although those forces 
give the subdominant contributions, some of the subleading terms in Eq. f)3.32p show resonant 
behaviors and become comparable to the leading term at u± = 0, U2 = or uj\ = U2. This is 
because the corresponding forces in the equation of motion Eq. A3. 12[) do not oscillate for these 
values of frequencies, that is, the forces are not averaged in the case of vortices with synchronized 
orientations. In the case of m e > m g with qi ■ q 2 — 0, the leading term Eq. f !3.34p vanishes and 
the terms with the resonant behavior become dominant as shown in Fig. [71 
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(a) a = 5 (b) a = 7 (c) a = 9 

Fig. 7: Scattering angles as a function of 72 for m g = m e = 1, o = 5, 7,9, 71 = 3, qi = n 2 o = (1,0,0) 
and q 2 = nio = (0, 1,0). The numerical results and those of free-motion approximation are denoted by 
dotted lines and solid lines, respectively. The resonances can be seen at 72 = and 72 = 71 = 3. 

3.4.2 Exchange of Q-charges 

In the scattering process, the charge of individual vortex in Eq. (13. 6 p is not conserved and its 
time dependence is given by the equation of motion for the orientation. Substituting the free 
motion Eq. (I3.27P into the equation of motion Eq. (I3.13P and integrating the right-hand side, we 
obtain the increment of the charge of the first vortex AQi (= —AQ 2 ) as 

AQi = ° 3 ™ gV (7172 qi x q 2 e~ m * a + 72 q\ x n 10 e ~ m ^V^ - 7l ft x n 20 e~ m ^ 

2n 2 c 2 m a v 1 + 2o , t79 rrn — ; — ^2 
^ ^7172 Rp -> 5 2 mgav / 1+(7l+72) 2 

9 2 v / l + (7i+72) 2 1 1 
2Ti 2 c 2 m a v l — 2^1 ^9 rr~/ vi 

9 2 V 1 + (7i " 7 2 ) 2 
Here, we have used the following integration formula 

vt ± ia 



-<-> 



(3.36) 



y/(vt) 2 + a 2 



dt( Ti = = ) K 2 {m^{vt) 2 + a 2 )e^ 
2tt 



(™) 2 ^ v/(mt0 2 + y 

For generic values of the angular velocities of the orientations 7/, the leading contribution is given 
by 

AQi « ^ 7i72 9i x <f 2 e""* . (3.38) 

g 2 

This contribution comes from the electric coupling Eq. ( 13. 18ft as in the case of Az\ 2 . The resonant 
behavior of the sub leading terms can also be seen at 71 = 0, 72 = and 71 = ±72. In particular, 
if the vortices are initially uncharged (71 = 72 = 0), the increment AQi is given by 

7r 2 c 2 m„v 

AQi = n 10 x n 20 e~ m * a . (3.39) 

g 2 

Therefore, the vortices are charged after the scattering even if Qj = initially. 
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3.4.3 Exchange of energy 



Let us now show that the energy can be transferee! between two non-Abelian vortices through the 
scattering. This is another characteristic feature of the non-Abelian vortices. Since two Abelian 
vortices are identical objects without any internal degree of freedom, the energy transfer never 
occurs in the scattering of two Abelian vortices since only the elastic scattering is possible for 
two identical vortices. 

The energy of an isolated vortex is given by the sum of the kinetic energies of the spatial 
motion and the orientation 

7TT77 2 

E I = ^\z I \ 2 + ^\Q I \ 2 . (3.40) 
g z Att 

Assuming that Ai 12 and AQi are sufficiently small, we can calculate the total change in energy, 
AE 1 (= -AE 2 ), as 



ixm 



2 „i 



AE 1 « —^(AiiZi + irjAlO + f-AQi-Qi (3.41) 
g 2 2vr 



g a v /T+7| 



-2nE c 2 glll2 (ft x q 2 ) ■ (n w e^v 7 ^? + n 20 e 
-2nE .2 (Ti -T2)(l + 2Ti7 ^ ) ^ . ^ e - mg a^M^f 

y/1 + (7l + 7 2 ) 2 



9 v / l + (7i-72) 2 



he energy is exchanged between 
3, while the resonances dominate 



where Eq = nm 2 g v 2 / (2g 2 ) is the kinetic energy associated to the initial relative velocity. Note 
that Z\ = \z\i in the center of mass frame. This shows that 
two vortices only through the subleading interactions in genera 
when 71 = 72 and 71 = (72 = 0), as shown in Fig. El 

3.5 Zero impact parameter 

In this section, we discuss the vortex scattering with zero impact parameter. For simplicity, we 
consider the effective Lagrangian restricted to the subspace of the moduli space given by 

z 12 = reR, /3 7 = e iw . (3.42) 

This subspace is a fixed point set of the following reflection symmetry on the moduli space 

zn -> 2i2, Pi -> TP- (3.43) 

Pi 



7 The leading terms are of order max{e meQ , e m s a }. 
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Fig. 8: The energy transfer AE\/Eq as a function of 72 (m g = m e = 1). (a) a = 9, 71 = 2, q± = q 2 = 
(1,0,0), n w = (0,1,0) and n 20 = (0,0,1), (b) a = 9, 71 = 1, n 10 = n 20 = (0,0,1), ft = (1,0,0) and 
<?2 = (0, 1, 0). The resonances can be seen at 72 = 71 in (a) and at 72 = in (b). 

According to the principle of symmetric criticality, any solution of the equation of motion re- 
stricted on such an invariant submanifold is automatically a stationary point of the original 
effective action. On this submanifold, the Lagrangian becomes 

L = —2 [l - 2c 2 e K (m e r) - 2c 2 g K (m g r) cos<p r ] m 2 r 2 

+ + ^g K o( m 9 r )0- + 3cosy?r)] y? 2 - ^c 2 g Ki(m g r) sin <p r m g fcp r 

? 

where y?o = + V2) an d = ¥>i — ^2- 111 this case, there exists only one non-vanishing 

conserved charge 

Itt 

which corresponds to the third component (Q)3 of the conserved Noether charge. By a Legendre 
transformation, tp can be eliminated from the effective Lagrangian in favor of Q. As a result, the 
following potential is induced by the Noether charge 

V Q = fV (l - 2^oM sin 2 I)"' « |V (l + 2c 2 g K (m g r) sin 2 |) . (3.46) 

This potential is shown in Fig. [9j Then, the equations of motion are given by 

■ 2 



+ 4 (l - 2c 2 g K (m g r) sin 2 ^) <pl (3.44) 



Q = ^ (l - 2c 2 g K (m g r) sin 2 ^) 00, (3.45) 



\c 2 jn e Ki{m e r) + c 2 m g Ki(m g r) cos<^ r ] 



r 2 / 4 



2c g K (m g r) sin ip r fip r + — ifi(m ff r) sin — ( —Q - ip r ) , (3.47) 



</3 r = — c g m g K 2 (m g r) sin ip r r + c g m g K 1 (m g r)(l + 3 cos cp r ) fip r 

-^-K (m g r) sinyv (j^Q 2 - 30 2 ) . (3.48) 
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If the charge Q is sufficiently large compared to (f r ,r and the relative angle (p r takes a generic 
value, two vortices recoil due to the potential Vq, in which case we can continue to use the 
asymptotic metric. This phenomenon can be seen in the numerical calculation Fig.[4]-(a). On the 
other hand, when Q is small or ip r m (ni • n 2 ~ 1), the vortices can closely approach along the 
valley of the potential Vq, in which case we can trace the dynamics until the asymptotic metric 
becomes invalid. Note that the cases with Q « and <p r ~ const, correspond to the resonance 
at 7i = 72. 
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Fig. 9: The induced potential VqKq Q /8tt) with c e = c g = 1.708 and m e = m g = 1. 



4 BPS Dyonic vortices 



4.1 Dyonic vortices from a mass deformation 



In this section we discuss the dynamics of BPS dyonic vortices in the U (2) gauge theory. In order 
to obtain the dyonic vortices, we add the adjoint scalar field E = H°t a + T, a t a and the following 
kinetic and mass terms [21] to the original Lagrangian Eq. (12. ip 

£kinadj = ^Tr [P^E^E] , £ m ass = Tr [(EH - HM)(EH - HMf] , (4.1) 

where M is a 2-by-2 mass matrix. This is the only consistent mass deformation preserving all 
supersymmetry, when embedded into the supersymmetric theory with eight supercharges^. The 
mass matrix can be written as a linear combination of the Pauli matrices 



M = —m ■ a. 
2 



(4.2) 



8 Although the adjoint scalar E and the mass M can be triplets of SU{2) C SO(A) r in (2+l)-dimensional 
A/" = 4 supersymmetric gauge theories, only one component of each triplets is relevant to the dyonic vortices. 
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In the following, we choose the mass to be m = (0, 0,m). This mass term induces the non- 
vanishing VEV for the adjoint scalar 

(E) = M, (4.3) 

and breaks the SU(2)c+f symmetry to U(l) (U(l) m from now on), for which the corresponding 
conserved charge is the component of Q parallel to rh. 

In this mass deformed model, the BPS bound for the energy density is given by [23| I4T] 



£ > ~c fu - ie ij diTr [VjHH*] + iTr [HMV H^ - V HMH^} + -^Tr [F 0l E] . (4.4) 

The first two terms are the energy density of vortices and the third term is the conserved Noether 
charge for the unbroken U(l) m symmetry. The forth term is the divergence of the the electric 
flux in the internal direction specified by the adjoint scalar E. We call this quantity simply 
"electric charge density" . The BPS bound is saturated if Eq. ( 12. 9 j) and the following equations 
are satisfied 

VqH — —i(EH — HM), T> i H = F 0i , T? E = 0, (4.5) 



2 Q 2 

e l g z 



-i 



(HVqH^ — VqHH^). (4.6) 



These BPS equations for dyonic vortices can be simplified by choosing the gauge in which the 
time component of the gauge field takes the form 

Wq = M - E. (4.7) 

Then, Eqs. (EOj]) . fET~6]) become 

doH = i[H,M\, d Wi = i[W h Ml SbE = i[E,M], (4.8) 



V{D, 



e l g 2 



{HH\E} = -2HMH ] . (4.9) 



From these equations we find that the BPS dyonic vortex solutions can be obtained by solving 
Eq. (14. 9 p with respect to E in a static vortex background satisfying Eq. (I2.9p and then rotating 
the orientation as 

H ->■ WHU, Wi ->■ rfWiU, E -> rfEU, U = e iMt e U(l) m . (4.10) 

Therefore, the dyonic vortex solutions are stable stationary configurations with rotating orien- 
tation. An important fact is that the no static force is exerted among the BPS dyonic vortices, 
that is, the electrostatic interaction Eq. ( 13 . 1 8[) is canceled by another interaction induced by the 
adjoint scalar E. 
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One can find the energy of the BPS saturated configurations by integrating the right-hand 
side of Eq. ( I4~4j) as 



E = 2-Kck + m ■ Q, 



where we have used 



m-Q = i J d 2 xTi [HMV H ] - V Q HMH ] ] . 



(4.11) 



(4.12) 



Note that the forth term in Eq. (I4.4p . as well as the second term, has no contribution to the 
total energy since the electric field is screened and decays exponentially in the Higgs phase. Note 
that the dyonic vortices have non-trivial charge distributions (see Fig. lTO]) . although they have 
no total electric charge. 
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(a) Tr[F 0i E] 



(b) ^Tr[F 04 S] 



(c) iTr[HMV H^ - (c.c.)] 



Fig. 10: An example of a single dyonic vortex configuration (numerical solution) for m g = 
m e = m = 1, 9q = | . (a) The electric flux radially spreads from the vortex center and decays 
exponentially, (b) The electric charge density is positive inside the vortex core while it is negative 
around the core, (c) The conserved Noether charge density is positive everywhere and the total 
charge is non-zero. 



We can also discuss the dyonic configurations in the effective theory of vortices. The defor- 
mation term Eq. 04. ip induces the following potential on the moduli space 



k l = m • C 



(4.13) 



where is the holomorphic Killing vector for SU(2) defined in Eq.( l3.5p and k % is that for U(l) m . 
In this deformed effective theory of vortices, the energy is bounded by the Noether charge Q 



E 



> ±m-Q. (4.14) 
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Therefore the solutions of the dyonic vortices, which are trajectories on the moduli space with 
minimum energy for a given value of m ■ Q, are determined from the BPS equation 

ft = ±k\ (4.15) 

In the rest of this paper, we take the positive sign for BPS dyonic vortices. In the case of 
m = (0, 0, m), the dyonic vortex solution is given by 

Zl = z I0 , ft = p I0 e mt , (V fc*A = im J2^A\ , (4.16) 

where zj and flio are complex constants which are related to the Noether charge by 

m-Q = g fj (ftk j + } j k i ) = 2g fj k i k j , (4.17) 

where we have used Eq. (I4.15p . Again, we find that the orientations of the vortices are rotating 
with the same period. 



Fig. 11: The trajectory of the internal orientation for a single dyonic vortex. The orientation 
rotates around a small circle at a fixed latitude 8q. 

As an example let us consider the case of a single dyonic vortex. Since the dynamics of the 
vortex position is trivial in this case, we consider only the vortex orientation. For a single vortex, 
the mass deformed effective Lagrangian can be written in terms of the spherical coordinates 

P = tan f e itp as 

6 2 + sm 2 6(if 2 -m 2 )] . (4.18) 



Q = — = —sm 2 6<j>- (4.19) 



The Noether charge Q = rh ■ Q/\m\ in this case is the conjugate momentum of <p 

dL 2tt 
g 

The dyonic vortex solution is given by 

2 TV 771 

e = 6 , <p = mt + <p , Q = — — sin 2 6 . (4.20) 

g 2 
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where 9q and tpo are constants. The trajectory corresponding to the dyonic vortex is a small 
circle at 9q, on which the following effective potential is minimized 

V{6) = V Q + V m = -^- Q Q 2 + ^P-m 2 > ^-sin 2 9 , (4.21) 

An sm 9 g z g z 

where Vq is the potential induced by the Noether charge Q. Let us consider a small fluctuation 

59 from the trajectory ( 14.20P corresponding to the dyonic vortex: 9 = 9q + 59. The effective 

potential for 59 is given by 

V(9) « sin 2 9 + ^- cos 2 9 59 2 . (4.22) 

9 9 

Therefore, if the dyonic vortex is excited, the angle parameter 9 oscillates with frequency u = 
2m cos Oq, which is generically of order m = \rh\. Thus we find that there exist massive modes 
around the dyonic vortex configuration. 



4.2 Dynamics of dyonic vortices 

In this section we discuss the dynamics of two dyonic vortices by using the equations of motion 
Eqs. (13.121) and (13.131) modified by the mass term. Since the action of the Killing vector on the 
moduli parameters is given by 

(m • £ + rh • £*) zi = 0, (m • £ + rh ■ £*) nj = rh x nj, (4.23) 

the Killing potential V m can be obtained from the kinetic terms of the effective Lagrangian by 
dropping Zi and replacing ni as 

nj — > rh x nj. (4.24) 

The full mass deformed Lagrangian can be found in Appendix |A] Then the BPS equations for 
dyonic vortices can be rewritten as 



n j = rh x nj. (4.25) 



The BPS solution takes the form 

1 



ni = — 
m l 



yjl — \vi\ 2 rh + cos(mt) mvj + sin(mt) rh x iJ[ , (4.26) 
where i>V are vectors such that iJj ■ rh = 0. 



Let us next see how the equations of motion for the relative position and the orientations 
are modified. The mass deformation changes d\ ■ a<i in the last term of the equation of motion 
Eq. (l3~T2l) as 

ol\ • c?2 - > c?i • c?2 — c?i • a 2 ) (4.27) 
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where a! l is the vectors obtained from aj = Hi — %Hi x Hi by replacing Hi with m x Hi 

a'i = m x Hi — iHi x (m x Hi). (4.28) 

Therefore, if the pair of the vortices are near BPS Hi ~ m x Hi, the contributions from the last 
term of Eq. f 1 3 . 1 2 j) and the mass deformation are small 

«i ■ a 2 - a[ ■ a' 2 « 0. (4.29) 

Similarly, we can also show that there exists a similar cancellation in the equation of motion for 
the orientation Eq. ( 13. 131) 

a.\{H2 ■ Si + 2H\ ■ 612) — «'i(^2 • ct\ + 2ni • q.' 2 ) ~ 0. (4.30) 

In the following, we will see that due to these cancellation, the behavior of the dyonic vortices 
are quite different from that of vortices before the mass deformation. 



4.2.1 Angular momentum, Lorentz Force and a Bound State 

Let us first discuss the case in which the relative velocity of the dyonic vortices are small. By 
using the effective Lagrangian Eq. (1A.1I) . we can write down the conserved angular momentum 
in terms of the radial coordinates 212 = re %x as 

dL Timl „ 

Lx = a£ « -fr*x + A x {r), (4.31) 

where we have neglected a small term proportional to x^o- The function A x is given by 

A x (r) = -(h 1 -Q 2 + H 2 - Qi) c\K x {m g r)m g r. (4.32) 

The first term in Eq. (14.311) is the ordinary angular momentum of a free particle and the additional 
term A x (r) can be interpreted as a "gauge potential" generated by the motion of the orientational 
moduli. 

In the massless case discussed in section [31 the contribution of A x (r) is small compared with 
the leading terms. This can be seen by averaging over the rapid motions of the orientations. 
Since they rotate around great circles of CP 1 , it follows that 

^H 1 -Q 2 + H 2 -Qi^ = 0. (4.33) 

The mass deformation drastically changes this situation. Since the orientations rotate around 
small circles in the case of the dyonic vortices, the average of the gauge potential is non-zero 

(nt ■ Q 2 + H 2 ■ Qi} = -j rh ■ (Hx + n 2 )(l - Hx ■ H 2 ). (4.34) 
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Note that m ■ ni and n\ ■ fi2 are independent of time for BPS dyonic vortices and slowly vary 
for an interacting pair of vortices. Due to the contribution from the gauge potential A x , the 
motions of the dyonic vortices with small relative velocity v = \z 12 \ becomes quite different from 
that without the mass deformation. For a pair of near BPS dyonic vortices nj m m x fij, the 
non-trivial gauge potential A x (r) gives the dominant contribution via the Lorentz force 

,2 



7rm: 



-Z12 



-iBz 12 , 



B(r) 



ldA x (r) 



r dr 

Assuming that the relative velocity v = \zn\ and the massive oscillations [59 in Eq. f)4.22p 
small, we can write down the energy of this system as 



(4.35) 



arc 



E 



irm 



2g 2 



9f2 



7l 



Mr) 



m g r 



(4.36) 



where we have neglected irrelevant terms. For an arbitrary relative distance r 
always adjust L x so that Vl = 0, namely 

L x = A x (r ). 



ro, we can 



(4.37) 



This minimum energy configuration corresponds to the BPS solution, for which the velocity of 
the relative position is zero r — x — 0- As shown in Fig. [12]- (a), this is a stable point of the 
positive semi-definite potential Vl, and thus the relative distance oscillates around the minimum 
for a sufficiently small excitation energy. The time dependence of the relative angle x can be 
determined from Eq. (I4.3ip . As a result, we can find that the dyonic vortices drift along a contour 





(a) (b) 

Fig. 12: (a) The potential Vl(t) representing the Lorentz force for m g ro = 4. (b) Examples 
of orbits (numerical solutions) for m g v/\m\ = 0.02 (solid line), 0.04 (dotted line), 0.16 (dashed 
line). Solid line shows an extended coil winding circularly around the origin. 



line of the magnetic field similarly to a charged particle in a background magnetic field. The 
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orbit of each dyonic vortex in the z-plane takes the form of an extended coil winding around a 
contour line (see Fig. [TBI (a)). In this sense, a pair of dyonic vortices can form a "bound state" 
with a large relative distance r ^> m~ 1 . The frequency of the oscillation around the stable 
point r = r is given by 

,2 



9 



ixm; 



■B(r ) oc K (m g r Q ) 



(4.38) 



and is generically quite small <C \m\ for well-separated vortices with m 9 r ^> 1. The radius 
of the coil r co ii, the velocity of the drift vd and the period of the large circular motion Tb are 
respectively estimated as 

4:n\u) h \r 



'"coil 



VTLgV 



T, 



13 



m a v* 



(4.39) 



|w L | 2|w L | 

where we have assumed that the relative velocity v = \z±2\ is sufficiently small. Note that the 
dyonic vortices with v larger than a certain critical value run away to infinity as illustrated in 
Fig.EKb). 



y " 





2000 4000 6000 8000 10000 



(a) Orbit of an extended coil around a circle (b) Deviation of energy from the BPS mass 

Fig. 13: An example of the bound state (numerical solution) for m g v/\rh\ = 0.016. (a) The orbit in 
z-plane. (b) The deviation of energy AE\ (divided by the BPS mass) as a function of time. 



In this analysis, we have neglected subleading terms in the equation of motion. This is valid 
if the excitation energy is sufficiently small since the subleading terms just change the potential 
VL(r) slightly for a small relative velocity v. We have also assumed that the motions of the 
orientations are kept near BPS hi m m x Hi and fluctuations of the massive modes S8i are small. 
These assumptions are justified if we take the near BPS initial condition with sufficiently small 
v, since the geodesic force which is independent of iu is almost cancelled by the potential term 
induced by the mass deformation. Therefore, we expect that the bound state is stable as long as 
the deviation of the energy of each vortex is sufficiently small. Fig.fTBTfb) shows an example of 
the deviation of the energy from the BPS mass during the period Tb. 
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Now let us discuss the origin of the Lorentz force Eq. ( 14. 35ft . Since dyonic vortices have both 
magnetic and electric fluxes, their asymptotic interaction can be well described by regarding 
them as point-like sources with both magnetic and electric moments. The static electric-electric 
and magnetic-magnetic interactions are canceled by the other scalar fields, so that the dominant 
interaction for a small relative velocity is the electric-magnetic interaction which is proportional 
to Zi2- Therefore, each slowly moving dyonic vortex feels the Lorentz force of the form Eq. (14.351) 
from the fluxes of the other dyonic vortex. 



4.2.2 Scattering of dyonic vortices 



In this section, we discuss the scattering of dyonic vortices by using the free motion approximation 
as in section [3j For simplicity, we restrict the initial conditions for the vortices to the BPS states, 
that is 

1 



Z12 = vt + ia, 



m = — 

m l 



yl — |t//| 2 m + cos(mi) Vj + sin(mi) m x Vj 



(4.40) 



Substituting this into the equation of motion Eq. (I3.12p modified by the mass deformation 
Eq. (I4.27p . we obtain the following equation for the relative position z\2 

vt — ia 



Zl2 



1 1 

-m e v c 



2 + a 2 



i^(m e vW + a 2 ) 



.2 2 



vt — ia 



—m„v c a 

y 9 /r..j.\ 2 



-.Krimg^^vtf + a 2 ) rii ■ n 2 



-2vic 2 g K (m g >y (vt) 2 + a 2 ) m ■ (ni + n 2 )(l - n x ■ n 2 ) 



(4.41) 



Here n\ ■ n2 and rh ■ ni are constant since we have assumed that the orientation modes satisfy 
the BPS equations ?ij = mx fij. Note that the forces vanish in the limit v — > due to the BPS 
properties of the dyonic vortices. The total change in the relative velocity Az±2 is 

Az 12 



2niv 



-c 2 e e-^ a + -cy^n, ■ U2 - c 2 e -^ >-^ + ^) ) 
2 2 y y 



rrigV 



(4.42) 



Similarly, we can calculate the total change in the Q-charge 



-A(m ■ Q 2 ) 



Tic g mgV 



(Hi x r?2) • fh 



2m ■ (ni + 712) 



rrigV 



(4.43) 



From these results, we find that Azu and A(rh-Qi) do not vanish even for small relative velocitj^] 
v « 0. This is because the Lorentz force (the last term in Eq. (14.4ip ) is proportional to v and its 
integrated effect is independent of v. Before the mass deformation, the differences Azu and AQi 



9 Although the free motion approximation would not be a good approximation for v ps 0, the results Eqs. (|4.42[) 
and (|4.43[) are qualitatively correct if the impact parameter a is sufficiently large. 



32 



vanish for a pair of slowly moving BPS vortices (v — > 0, 7/ = 0) as we can see from Eqs. (13.331) 
and (I3.36p . Note that in both cases the interaction vanish in the limit v — > because of the BPS 
properties. This is one of the typical differences between BPS vortices before and after the mass 
deformation. 

Since the first two terms in Eq. (14.421) can be obtained by taking the limit 7^2 — > in 
Eq. (13.32 p . they consistently reduce to Eq. (I3.32p in the limit \m\ — > 0. Note that 77 oc \Qi\ are of 
order \m\ for the BPS dyonic vortices. The effect of the mass deformation lies in the last term, 
which becomes dominant if m g v \m\. On the other hand, the free-motion approximation gives 
quantitatively precise results when the relative velocity v is sufficiently large : m 9 v » me^. 
In such situation^, we can show that the impact parameter a should be sufficiently larger than 
the minimum of the potential tq which is given by 

KM = |L X | = ^(m g a)(m g v). {AAA) 

Therefore Eqs. (14.421) and (14.431) are valid only for the scattering of dyonic vortices repelled by 
the outer barrier of the potential V^ir). 

As we can see from Eqs. (I4.42p and (14.431) . there is no contribution of order m 2 oc (7/) 2 . 
This is one of the most striking features of the dyonic vortices, which makes the Lorentz force 
dominant and ensures the stability of the bound states. Strictly speaking, we have to consider 
contributions from the massive oscillations which we have ignored for simplicity. We expect that 
massive modes are irrelevant for the near BPS configurations, while they would cause instabilities 
of highly excited bound statei"! 



5 Conclusion and Discussion 

We have studied the dynamics of non-Abelian vortices by the moduli space approximation using 
the asymptotic metric of well-separated non-Abelian vortices in the U(N) gauge theory with 
N Higgs scalar fields in the fundamental representation. Since non-Abelian vortices carry the 
orientational moduli CP^ -1 individually, they can have Noether charges as internal momentum 
of them, whereas ANO vortices in the Abelian-Higgs model have no such charges. 

We have found that the vortices with the same charges repel while those with the opposite 
charges attract. We have shown that the charges of vortices can change during the scattering pro- 
cess with the total energy conserved; the kinetic energy of orientational moduli can be transfered 

10 In the opposite case m g v 3> |m|, the approximation is always valid, but the effects of the mass deformation 
become subleading. For instance, there is no minimum in the potential Vh(r) in such a case. 

11 In the latter case, there would be corrections with resonant behaviors in the scattering angle Ai 12 , etc. 
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to the kinetic energy in real space and vice versa. As results, we have found that in scattering of 
two non-Abelian vortices, i) the scattering angle depends on the internal orientation, especially 
parallel orientations give repulsion while anti-parallel orientations give attraction, ii) the energy 
of real and internal spaces can be transfered, iii) the energy and charge transfer between two 
vortices occur, and iv) some resonances appears due to synchronization of the orientations. 

By introducing the mass deformation into the original theory, the color-flavor symmetry 
SU(N) is explicitly broken to [/(l)^ -1 . Noether charges of these [/(l) 1 ^' 1 can be actually regard 
as [/(l)^ -1 gauge charges where the Cartan subgroup of U(N) survives in the low energy theory 
in the large mass limit. We have shown that the dominance of the Lorentz force between vortices 
gives a bound state with a coiling orbit. 

Here we address several discussions. 
One may question what are the conditions for the existence of bound states in more general 
cases. To form a bound state, therefore, vortices are needed only to have 

1. a magnetic flux (a topological charge) as the definition, 

2. a sufficiently large conserved charge, 

3. a sufficiently small static force. 

Since large conserved charges generically induce large interactions between vortices, we need to 
prepare a special system for cancellation of static forces like a BPS state. Even in the Abelian 
case, vortices can have conserved charges. For instance, Abelian semi-local vortices [15*1 HB] can 
acquire Q-charges after the mass deformation. In the strong coupling limit e — > oo, they reduce 
to Q-lumps and their interactions have been studied [T9| Wf\. 

As future problems, the followings might be interesting. 

Our study in this paper has been based on the moduli space metric for well-separated vortices 
and therefore we cannot discuss vortices with a small separation. One question is what the fate 
of a pair of vortices with anti-parallel charges is. Will they form a bound state with a very small 
separation? Or will they reduce to a composite state of coincident vortices (singlet or triplet) 
[35l [38] after radiations of massive particles? The head-on collision of two vortices was studied 
previously based on the metric on the moduli subspace of two coincident vortices (at zero dis- 
tance) [39] . However the moduli space metric for two vortices at arbitrary distance is not known, 
which is needed to answer the above questions. Also, beyond the moduli space approximation, 
one has to study dynamics by numerically solving the original equations of motion. 
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We have found that energy transfer and Noether charge transfer between two scattering 
vortices. This property may support the Boltzman's principle of equality when one considers 
statistical mechanics of many vortices. Partition function of a gas of non-Abelian vortices was 
calculated in [3S| with assuming the principle of equality so that the calculation is reduced to 
the integration over the moduli space of vortices. Therefore, one important question is if a 
large number of scattering of vortices leads to the ergodic theorem so that one can assume the 
Boltzman's principle of equality. 

Non-Abelian vortices are called semi-local if there are more flavors than the number of color. 
Non-Abelian semi-local vortices have size moduli and they reduces to local vortices when the 
size modulus is sent to zero. One of characteristic properties of semi-local vortices is that their 
profile functions decay polynomially but not exponentially as local vortices. Consequently the 
size modulus is non-normalizable. Moreover the CP N ~ 1 orientational moduli of a single non- 
Abelian semi-local vortex are also non-normalizable [19] when the size modulus is non-zero, but 
they are normalizable only when the size modulus is zero [SD]. No metric can be defined for 
non-normalizable moduli, and dynamics cannot be discussed for those. However we can discuss 
the dynamics of two semi-local vortices because relative orientations and relative size between 
them are normalizable. In fact dynamics such as head-on collision was studied for Abelian 
[IB] and non-Abelian [3H] semi-local vortices. Since semi-local vortices can be approximated by 
lump solutions at large distance compared with Compton wave length of massive particles, the 
dynamics of semi-local vortices at large distance can be approximated by that of lumps |19j . 

One interesting generalization is changing geometry from a flat space to geometry with non- 
trivial cycles such as a cylinder [51] , a torus [H] , and Riemann surfaces with higher genus [55] . For 
semi-local vortices in compact Riemann surfaces, there is no problem of the non-normalizability 
and we can discuss their dynamics without introducing any cut-off scale by hand. In fact the 
moduli space metric was found for well-separated non-Abelian vortices on Riemann surfaces 
[56] • However dynamics during a long time will be difficult using the asymptotic metric because 
distances between vortices cannot be kept large for compact spaces. 

Another interesting extension is changing gauge groups. Non-Abelian vortices were extended 
to arbitrary gauge groups G in the form of (U(l) x G)/C(G) with the center C(G) of G [5Tj . 
Especially the cases of G = SO(N), USp(2N) have been studied in detail [52], [53] . However they 
are semi-local vortices in general [52] so that dynamics at large separation is the one of lumps 

P2I- 

We have studied dynamics of BPS vortices. There exist static force between non-BPS vortices 
[57] . Superconductors are classified into two types, type I and type II, depending on whether 
the static force between two vortices is attractive(type I) or repulsive (type II). A classification 
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of non-Abelian superconductors will be more complicated due to the existence of charges of non- 
BPS non-Abelian vortices [58]. When the couplings are close to the critical coupling (near BPS), 
dynamics of non-BPS vortices can be studied by the moduli space approximation plus a potential 
for the static forces, which remains as an important future problem. Inclusion of Chern-Simons 
terms [59] is also possible extension. 

Our studies have been restricted to particle dynamics, namely vortices in 2+1 dimensions. In 
3+1 dimensions, vortices are (cosmic) strings which have one spatial dimension in their world- 
sheet. The moduli space approximation can be applied when the angle between two cosmic 
strings is small. Collision of two non-Abelian cosmic strings was studied based on the metric 
on the moduli subspace of two coincident vortices [501 133- Especially it was found in [39J that 
orientational moduli of two non-Abelian vortices must be aligned and scatter with 90 degree 
angle, when they collide in head on, except for a fine-tuned collision. This implies that two 
non-Abelian cosmic strings reconnect each other. This result was obtained just before and after 
the collision moment in the linear order in time. On the other hand, Eq. (I3.26P shows that 
orientational moduli of two vortices in head-on collision tend to be aligned. It suggests that two 
non-Abelian cosmic strings consistently reconnect each other as a long time behavior by feeling 
attraction between two orientational moduli like a ferromagnet. However we also have found in 
this paper that two vortices repel when charges induced by the motion of internal orientational 
moduli of two vortices are opposite. In this case, we expect that two cosmic strings merely 
reconnect each other. 

We have found many new feature of dynamics or scattering of non-Abelian vortices, which are 
based on the existence of non-Abelian internal moduli of individual solitons. Therefore the similar 
properties should hold for scattering of other kinds of non-Abelian solitons such as Yang-Mills 
instantons, non-Abelian monopoles [5U], and non-Abelian kinks |32j . 
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A The action in terms of the orientation vectors Hj 



In the case of U(2) gauge theory, it is convenient to describe orientations of vortices in terms of 
three-component unit vectors Hi. The Kahler potential in Eqs. ( 12. 17ft and (12. 18ft give the effective 
Lagrangian 

L = Lkin + L int — V m . (A.l) 

Introducing Lagrange multipliers Xj G R, the effective Lagrangian C c s can be rewritten in terms 
of Hi. The kinetic terms giving free motions are obtained as 



E 7 (\ m2 M + l&l 2 + MN 2 - !)) ■ ( A - 2 ) 



Note that the quantity G12 in the Kahler potential is G12 — H\-H 2 and 

9 L 1 H \ 9 

d/3— — - ( dH — tz^(H ■ dn) — in x dn) ■ — . (A. 3) 

op 2 \ \n\ 2 J On 

Therefore kinetic terms L- mt describing interactions between vortices can be directly calculated 
from the Kahler potential as 

g 2 1 

— L int = _ 2 m 3 V + { c l K o( m e r ) + • H 2 )K (m g r)) 

+c 2 I ((Ht ■ H 2 ) + (H 2 ■ Hi) J rrigf + (Hi x H 2 ) ■ {n x - n 2 )m 9 rx j Kx(m g r) 

-c 2 i—Hi ■ n 2 (\ni\ 2 + |n 2 | 2 ) + n x ■ n 2 + (H-y x H x ) ■ (H 2 x H 2 )\ K (m g r), (A.4) 

where z\ — z 2 = z\ 2 = re %x . The potential V m induced by the mass deformation can be obtained 
by just replacing zi — > and Hi — > m x Hi as 

2^ m = ^ x ^j| 2 ~ c q| ~ H'i " H 2 (\rh x ni| 2 + \rh x n 2 \ 2 ) 

+ (m x r?i) • (m x n 2 ) + [H\ x (m x ni)) ■ (n 2 x (m x n 2 )) j-L'o(wigr). (A. 5) 

Due to the Lagrange multipliers, the equations of motion for the orientations become 

M _ ^(9C d dC\ n _ fdC d dC\ n . . r . 

^-^{mrJtWj^- * n,x {^-7tW,) =0 - (A ' 6) 

For instance, the equation of motion of the orientation n for a single vortex is given by 

H x H + (m • H) fa x H = 0, (A. 7) 

with the conservation conditions H ■ H = 0, |n| 2 = 1. 
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